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Abstract 

In this paper, we complement recent results of Bronski and Johnson and of John- 
son and Zumbrun concerning the modulational stability of spatially periodic traveling 
wave solutions of the generalized Korteweg-de Vries equation. In this previous work it 
was shown by rigorous Evans function calculations that the formal slow modulation ap- 
proximation resulting in the Whitham system accurately describes the spectral stability 
to long wavelength perturbations. Here, we reproduce this result without reference to 
the Evans function by using direct Bloch-expansion methods and spectral perturbation 
analysis. This approach has the advantage of applying also in the more general multi- 
periodic setting where no conveniently computable Evans function is yet devised. In 
particular, we complement the picture of modulational stability described by Bronski 
and Johnson by analyzing the projectors onto the total eigenspace bifurcating from the 
origin in a neighborhood of the origin and zero Floquet parameter. We show the result- 
ing linear system is equivalent, to leading order and up to conjugation, to the Whitham 
system and that, consequently, the characteristic polynomial of this system agrees (to 
leading order) with the linearized dispersion relation derived through Evans function 
calculation. 



1 Introduction 



In this paper, we consider traveling wave solutions of the Korteweg-de Vries equation 

(1-1) U t = Uxxx + f(u)x 
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where u is a scalar, x,t £ R, and / G C 2 (R) is a suitable nonlinearity. Such equations 
arise in a variety of applications. For example, in the case f(u) = u 2 equation (jl.ip 
reduces to the well known KdV equation, which serves as a canonical model of unidirectional 
weakly dispersive nonlinear wave propagation [KdV] ISH] ■ Furthermore, the cases f(u) = 
/3-u 3 , /3 = ±1, correspond to the modified KdV equation which arises as a model for large 
amplitude internal waves in a density stratified medium, as well as a model for Fermi- 
Pasta-Ulam lattices with bistable nonlinearity |BM| |BOj . In each of these two cases, the 
corresponding PDE can be realized as a compatibility condition for a particular Lax pair 
and hence the corresponding Cauchy problem can (in principle) be completely solved via 
the famous inverse scattering transform^. However, there are a variety of applications 
in which equations of form (jl.lj) arise which are not completely integrable and hence the 
inverse scattering transform can not be applied. For example, in plasma physics equations of 
the form (jl.ip arise with a wide variety of power-law nonlinearities depending on particular 
physical considerations [KD] [M] [MS]. Thus, in order to accommodate as many applications 
as possible, the methods employed in this paper will not rely on complete integrability of 
the PDE (jl.ip . Instead, we will make use of the integrability of the ordinary differential 
equation governing the traveling wave profiles. This ODE is always Hamiltonian, regardless 
of the integrability of the corresponding PDE. 

It is well known that (jl.ip admits traveling wave solutions of the form 

(1.2) u(x, t) = u c (x + ct) 

for c > 0. Historically, there has been much interest in the stability of traveling waves of 
the form (1.2) where the profile u c decays exponentially to zero as its argument becomes 
unbounded. Such waves were originally discovered by Scott Russell in the case of the 
KdV where the traveling wave is termed a soliton. While (jl.ip does not in general possess 
soliton solutions, which requires complete integrability of the governing PDE, it nonetheless 
admits exponentially decaying traveling wave profiles known as solitary waves. Moreover, 
the stability of such solitary waves is well understood and dates back to the pioneering 
work of Benjamin [B], which was then further developed by Bona, Frillakis, Grillakis et. 
al., Bona et. al, Pego and Weinstein, Weinstein, and many others. 

The subject of this paper, however, is the stability of the traveling wave solutions of 
(jl.ip which are periodic, i.e. solutions of the form (jl.2p such that u c is T-periodic. Unlike 
the corresponding solitary wave theory, the stability of such periodic traveling waves is 
much less understood. Existing results in general fall into two categories: spectral stability 
with respect to localized or periodic perturbations [Br Jj IBP] . or nonlinear (orbital) stability 
with respect to periodic perturbations [ BrJKj IT U IDKj . In the present manuscript, we will 
consider spectral stability to long-wavelength localized perturbations, i.e. spectral stability 
to slow modulations of the underlying wave. Such an analysis was recently carried out by 
Bronski and Johnson in [BrJj by Evans function techniques. In particular, by studying the 
Evans function in a neighborhood of the origin a modulational stability index was derived 

1 More precisely, the Cauchy problem for equation p. II) can be solved via the inverse scattering transform 
if and only if / is a cubic polynomial. 
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which uses geometric information concerning the underlying wave to determine the local 
structure of the I? spectrum of the linearized operator at the origin. This sign of this index 
was then evaluated in several examples in the subsequent work of Bronski, Johnson, and 
Kapitula [BrJKj . 

Our approach is essentially different from [Br J]. In particular, we study the modulational 
stability by using a direct Bloch-expansion of the linearized eigenvalue problem as opposed 
to more familiar Evans function techniques. This strategy offers several novel advantages. 
First off, the Evans function techniques of |BrJ] do not extend to multiple dimensions in 
a straight forward way. In particular, these techniques seem not to apply to study the 
stability of multiply periodic solutions, that is, solutions which are periodic in more than 
one linearly independent spatial directions. Such equations are prevalent in the context of 
viscous conservation laws in multiple dimensions. In contrast, the Bloch expansion approach 
of the current work generalizes to all dimensions with no difficulty. Secondly, the low- 
frequency Evans function analysis required to determine modulational instability is often 
quite tedious and difficult; see, for example, the analysis in the viscous conservation law case 
in |OZH |QZ3"1 IOZ44 ISelj . In contrast, however, we find that the Bloch-expansion methods of 
this paper are much more straight forward and reduces the problem to elementary matrix 
algebra. Using this approach, then, we reproduce the modulatioanl stability results of 
Bronski and Johnson with out any specific mention or use of the Evans funcition. While 
this direct approach has been greatly utilized in studying the stability of small amplitude 
solutions of (jl.ip (see |DK1 IGHl IHal IHK1 IHaLS| for example), we are unaware of any 
previous applications to arbitrary amplitude solutions. 

Given the above remarks, we choose to carry out our analysis first in the gKdV case since, 
quite surprisingly, all calculations can be done very explicitly and can be readily compared 
to the results of |BrJ] , Our hope is that the present analysis may serve as a blueprint of how 
to consider the stability multiply periodic structures in more general, and difficult, equations 
than simply the gKdV equation. We should remark, however, that a slight variation of our 
approach was attempted in the latter half of the analysis in [Br J]. However, the analysis 
was quite cumbersome requiring the asymptotic tracking of eigenfunctions which collapse 
to a Jordan block. Our approach does not require such delicate tracking and it is seen that 
we only need to asymptotically construct the projections onto the total eigenspace, which 
requires only tracking of an appropriate basis not necessarily consisting of eigenfunctions. 

Finally, our approach connects in a very interesting way to the recent work of Johnson 
and Zumbrun [JZl]. There, the authors studied the modulational stability of periodic gKdV 
waves in the context of Whitham theory; a well developed (formal) physical theory for 
dealing with such stability problems. The calculation proceeded by rescaling the governing 
PDE via the change of variables (x, t) i— > (ex, et) and then uses a WKB approximation of the 
solution to find a homogenized system which describes the mean behavior of the resulting 
approximation. In particular, it was found that a necessary condition for the stability of such 
solutions is the hyperbolicity i.e., local well-posedness of the resulting first order system 
of partial differential equations by demonstrating that the characteristic polynomial of the 
linearized Whitham system accurately describes, to first order, the linearized dispersion 
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relation arising in the Evans function analysis in [BrJj. As seen in the recent work of [JZ4] 
and [JZ5j . there is a deeper analogy between the low-frequency linearized dispersion relation 
and the Whitham averaged system at the structural level, suggesting a useful rescaling of the 
low-frequency perturbation problem. Moreover, it was seen that the Whitham modulation 
equations can provide invaluable information concerning not only the stability of nonlinear 
periodics, but also their long-time behavior. It is this intuition that motivates our low- 
frequency analysis, and ultimately leads to our derivation of the Whitham system for the 
gKdV using Bloch-wave expansions. In particular, this observation provides us with a 
rigorous method of justifying the Whitham modulation equations which is independent of 
Evans function techniques used in |JZ1| . As mentioned above, such a result is desirable 
when trying to study the stability of multiply periodic waves where one does not have a 
useful notion of an Evans function. 

Therefore, the goal of this paper is to reproduce many of the results of [BrJ| and [JZ1] 
in a method which is independent of the Evans function and can be applied to higher 
dimensional and more general settings with little extra effort. Of particular importance is 
the fact that we can rigorously justify the Whitham modulation, a construct which exists 
in any dimension, using more robust methods than the restrictive Evans function approach. 
The plan for the paper is as follows. In the next section we discuss the basic properties of the 
periodic traveling wave solutions of the gKdV equation (jl.ip . including a parametrization 
which we will find useful in our studies. We will then give a brief account of the results and 
methods of the papers |BrJ] and [JZlJ concerning the Evans function and Whitham theory 
approaches. We will then begin our analysis by discussing the Bloch decomposition of the 
linaerized operator about the periodic traveling wave. In particular, we will show that the 
projection of the resulting Bloch eigenvalue problem onto the total eigenspace is a three- 
by-three matrix which is equivalent (up to similarity) to the Whitham system. We will 
then end with some concluding remarks and a discussion on consequences / open problems 
inferred from these results. 

2 Preliminaries 

Throughout this paper, we are concerned with spatially periodic traveling waves of the 
gKdV equation (jl.ip . To begin, we recall the basic properties of such solutions: for more 
information, see |BrJ] or [Jl| . 

Traveling wave solutions of the gKdV equation with wave speed c > are stationary 
solutions (II. ip in the moving coordinate frame x + ct and whose profiles satisfy the traveling 
wave ordinary differential equation 

(2.1) u xxx + f(u) x - cu x = 0. 

Clearly, this equation is Hamiltonian and can be reduced to quadrature through two inte- 
grations to the nonlinear oscillator equation 



(2.2) 



^ = E + au + \u 2 - F(u) 
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where F = f and F(0) = 0, and a and E are constants of integration. Thus, the existence 
of periodic orbits of (|2.ip can be trivially verified through phase plane analysis: a necessary 



and sufficient condition is that the effective potential energy V(u; a, c) := F(u) — ^u 2 — au 
have a local minimum. It follows then that the traveling wave solutions of (jl.ip form, up 
to translation, a three parameter family of solutions which can be parameterized by the 
constants a, E, and c. In particular, on an open (not necessarily connected) subset T> 
of M 3 = {(a,E,c) : a, E, c 6 M} in which the solutions to (I2.ip are periodic with period 
T = T(a,E,c): the boundary of these sets corresponds to homoclinic/heteroclinic orbits 
(the solitary waves) and equilibrium solutions. 

Moreover, we make the standard assumption that there exists simple roots u± with 
U- < u + such that E = V(u±;a,c) and E > V(u;a,c) for all u G (u_,u + ). It follows 
then that the roots u± are C 1 functions of the traveling wave parameters a, E, and c and 
that, without loss of generality, u(0) = u_. Under these assumptions the functions u± are 
square root branch points of the function \JE — V(u; o, c). In particular, the period of the 
corresponding periodic solution of (|2.ip can be expressed through the integral formula 



(2.3) T = T(a,E,c) = — f ^= du 

K ' 1 ^ 2 J r y/E-V(u;a,c) 

where integration over T represents a complete integration from u- to u+ , and then back to 
U- again: notice however that the branch of the square root must be chosen appropriately in 
each direction. Alternatively, you could interpret T as a loop (Jordan curve) in the complex 
plane which encloses a bounded set containing both u- and u+. By a standard procedure, 
the above integral can be regularized at the square root branch points and hence represents 
a C 1 function of the traveling wave parameters on V. Similarly, the conserved quantities of 
the gKdV flow can be expressed as 



(2.4) M = M(a,E,c) 



u du 



V2 r _ 

2 Jr .JE -V(u;a,c) 



„ „/ „ n v2 f u 2 du 

(2.6) n-HMQ-g f E -^tf- FM 

^ Jr \/E — V(u; a, c) 

where again these quantities, representing mass, momentum, and Hamiltonian, respectively, 
are finite and C 1 functions on T>. As seen in |BrJ[ IBrJKl Ulj. the gradients of the functions 
T, M, P : T> — > M play an important role in the modulational stability analysis of periodic 
traveling waves of the gKdV equation (jl.ip . 

Remark 2.1. Notice in the derivation of the gKdV [BaMo], the solution u can represent 
either the horizontal velocity of a wave profile, or the density of the wave. Thus, the function 
M : D — > M can properly be interpreted as a "mass" since it is the integral of the density 
over space. Similarly, the function P : D — > M can be interpreted as a "momentum" since 
it is the integral of the density times velocity over space. 
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To assist with calculations involving gradients of the above conserved quantities consid- 
ered as functions on P C K 3 , we note the following useful identity. The classical action (in 
the sense of action-angle variables) for the profile equation (|2.ip is given by 



K = <j> u x du = y/2 (j> y/E — V(u; a, c)du 



where the contour T is defined as above. This provides a useful generating function for 
the conserved quantities of the gKdV flow restricted to the manifold of periodic traveling 
waves. Specifically, the classical action satisfies 

, , dK dK dK 

as well as the identity 

K = H + aM + -P + ET. 
2 

Together, these relationships imply the important relation 

EV a , E , C T + aV a , E ,cM + ^V a ,E,cP + Va,E,cH = 0. 

So long as E ^ then, gradients of the period, which is not itself conserved, can be 
interchanged with gradientes of the genuine conserved quantities of the gKdV flow. As 
a result, all gradients and geometric conditions involved in the results in this paper can 
be expressed completely in terms of the gradients of the conserved quantities of the gKdV 
flow, which seems to be desired from a physical point of view. However, as the quantities T, 
M, and P arise most naturally in the analysis, we shall state our results in terms of these 
quantities alone. 

We now discuss our main assumptions concerning the parametrization of the family of 
periodic traveling wave solutions of (jl.ip . To begin, we assume throughout this paper the 
period is not at a critical point in the energy, i.e. that Te ^ 0. In other words, we assume 
that the period serves as a good local coordinate for nearby waves on T> with fixed wave 
speed c > and parameter a. As seen in |BrJKl |JT] . such an assumption is natural from the 
viewpoint of nonlinear stability. Moreover, we assume the period and mass provide good 
local coordinates for the periodic traveling waves of fixed wave speed c > 0. More precisely, 
given (ao,Eo,co) E V with Co > 0, we assume the map 

(a,E) ' y (T(a,E,c ),M(a,E,c )) 

have a unique C 1 inverse in a neighborhood of {uq,Eo) G M 2 , which is clearly equivalent 
with the non- vanishing of the Jacobian 

o(a, E) 
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at the point (ao, Eq, cq). As such Jacobians will be prevalent throughout our analysis, for 
notational simplicity we introduce the following Poisson bracket notation for two-by-two 
Jacobians 



and the corresponding notation {f,g,h} Xt y jZ for three-by-three Jacobians. Finally, we as- 
sume that the period, mass, and momentum provide good local coordinates for nearby 
periodic traveling wave solutions of the gKdV. That is, given (ao,-Eb>co) G T> with cq > 0, 
we assume that the Jacobian {T, M, P} a ,E,c is non-zero. While these re-parametrization 
conditions may seem obscure, the non-vanishing of these Jacobians has been seen to be 
essential in both the spectral and non-linear stability analysis of periodic gKdV waves in 
[BrJl [JT | IBrJKj . In particular, these Jacobians have been computed in [BrJKj for several 
power-law nonlinearities and, in the cases considered, has been shown to be generically 
non-zero. Moreover, such a non-degeneracy condition should not be surprising as a similar 
condition is often enforced in the stability theory of solitary waves (see [Bo} |B| IPW| ) . 

Now, fix a (ao, i^Oj Co) G *D. Then the stability of the corresponding periodic traveling 
wave solution may be studied directly by linearizing the PDE 

(2.8) u t = u xxx + f(u) x - cu x 

about the stationary solution it(-, ao, Eq, cq) and studying the L 2 (M) spectrum of the asso- 
ciated linearized operator 

d x C[u] ~d x {-d 2 x -f'(u)+c). 

As the coefficients of C[u] are T-periodic, Floquet theory implies the L 2 spectrum is purely 
continuous and corresponds to the union of the L°° eigenvalues corresponding to considering 
the linearized operator with periodic boundary conditions v(x + T) = e lK v(x) for all 
where k £ [— ir, it] is referred to as the Floquet exponent and is uniquely defined mod 2ir. In 
particular, it follows that \x G spec(d x C[u]) if and only if d x C[u] has a bounded eigenfunction 
v satisfying v{x + T) = e lK v(x) for some k£1. More precisely, writing the spectral problem 
d x C[u]v = fiv as a first order system 

(2.9) Y'(x)=H(x,n)Y(x), 

and defining the monodromy map M(/u) := <J>(T) < I>(0)~ 1 , where $ is a matrix solution of 
(12. 9p . it follows that \i G C belongs to the L 2 spectrum of d x C[u] if and only if the periodic 
Evans function 

D(p,K) := det (M(/x) -e iR l) 

vanishes for some k£1. 

When studying the modulational stability of the stationary periodic solution u(-; ao,Eo, Co) 
it suffices to study the zero set of the Evans function at low frequencies, i.e. seek solutions 
of D(fi,K) = for |(/i, k)\ 1. Indeed, notice that the low-frequency expansion h(k) for 
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(fi, re) near (0, 0) may be expected to determine the long-time behavior and can be derived 
by the lowest order terms of the Evans function in a neighborhood of (0, 0). As a first step 
in expanding D, notice by translation invariance of (I2.8D it follows that 

d x C[u]u x = 0, 

that is, u x is in the right T-periodic kernel of the d x C[u\. It immediately follows that 
D(0, 0) = 0, and hence to determine modulational stability we must find all solutions of the 
form (//(re), re) of the equation D(/j,,k) = in a neighborhood of re = 0. Using Noether's 
theorem, or appropriately differentiating the integrated profile equation (|2.2p . it follows that 
the functions u a and ue formally satisfy 

d x C[u]u a = d x C[u]uE = 0, d x C[u]u c = -u x . 

However, these functions are not in general T-periodic due to the secular dependence of the 
period on the parameters a, E, and c. Nevertheless, one can take linear combinations of 
these functions to form another T-periodic null-direction and a T-periodic function in a Jor- 
dan chain above the translation direction. A tedious, but fairly straightforward, calculation 
(see [BrJj) now yields 

(2.10) D(n, «) = A(/i, re) + 0(|/i| 4 + re 4 ), 

where A represents a homogeneous degree three polynomial of the variables \x and re. Defin- 
ing the projective coordinate y = — in a neighborhood of a = 0, it follows the modulatioanl 

LA 

stability of the underlying periodic wave u(-; do, E , cq) may then be determined by the 
discriminant of the polynomial R(y) := /i~ 3 A (1, —iy), which takes the explicit form 

(2.11) R(y) = -y 3 + | ({T, P} E , C + 2{M, P} a , E ) - \{T, M, P} a , E , c . 

Modulational stability then corresponds with R having three real roots, while the presence 
of root with non-zero imaginary part implies modulational instability. 

On the other hand, we may also study the stability of the solution u{-\ ao,Eo, cq) through 
the formal approach of Whitham theory. Indeed, recalling the recent work of [JZ1] we find 
upon rescaling by (x, t) t— > (ex, et) and carrying out a formal WKB expansion as e — > 
a closed form system of three averaged, or homogenized, equations of the form 

(2.12) d t (Mw, Poo, u) (u) + d x (F, G, H) (u) = 

where oj = T" 1 and it represents a periodic traveling wave solution in the vicinity of the 
underlying (fixed) periodic wave «(•; ao, Eq, cq). The problem of stability of u(-\ do, To, Co) to 
long-wavelength perturbations may heuristically expected to be related to the linearization 
of (|2.12p about the constant solution ii = u(-; do, Eq, Co), provided the WKB approxima- 
tion is justifiable by stability considerations. This leads one to the consideration of the 
homogeneous degree three linearized dispersion relation 

x, , , / d(Mu,Puj,uj) ik d (F, G, H) \ . . n 

(2.13) A(/i, re) := det U y ' ' ; - - V ^' > ) («(•; a ,E , c )) = 
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where \x corresponds to the Laplace frequency and k the Floquet exponent. 

The main result of [JZ1] was to demonstrate a direct relationship between the above 
approaches. In particular, the following theorem was proved. 

Theorem 1 ( |JZlj ). Under the assumptions that the Whitham system (|2.12p is of evolu- 
tionary type, i.e. ^ M q^ ,w ^ is invertible at (oo, Eq,cq) G T>, and that the matrix q^e c) * s 
invertible at (ao, Eq, cq) £ T>, we have that in a neighborhood o/(/i, k) = (0, 0) the asymptotic 
expansion 

D(fi,K) = T A(^ K ) + O(\fi\ 4 + K 4 ) 

for some constant r ^ 0. 

Remark 2.2. The assumption in Theorem[l}that d<yM Q^^ is invertible forces the Whitham 
system (|2.12p to be of evolutionary type. Moreover, notice that the Whitham system is in- 
herently a relation of functions of the variable u, while the Evans function calculations from 
[Br J] utilize a parametrization of V by the parameters (a,E,c). In order to compare the 
two linearized dispersion relations A(/x, k) and A(fx,K) then we must ensure that we may 
freely interchange between these variables, i.e. we must assume that the matrix g^jg c) * s 
invertible at the underlying periodic wave. 

That is, up to a constant, the dispersion relation (I2.13P for the homogonized system 
(|2.12p accurately describes the low-frequency limit of the exact linearized dispersion relation 
described in (I2.10p , As a result, Theorem [T] may be regarded as a justification of the WKB 
expansion and the formal Whitham procedure as applied to the gKdV equation. The 
importance of this result stems from our discussion in the previous section: although the 
Evans function techniques of [BrJj do not extend in a straight forward way to multiperiodic 
waves, the formal Whitham procedure may still be carried out nonetheless. However, it 
follows that we must find a more robust method of study to justify the Whitham expansion 
in higher dimensional cases. The purpose of this paper is to present precisely such a method 
using Bloch-expansions of the linearized operator near zero-frequency. Indeed, we will show 
how this general method in the case of the gKdV equation may be used to easily rigorously 
reproduce the linearized dispersion relation A(/x, k) corresponding to the Whitham system 
as well as justifying the Whitham expansion beyond stability to the level of long-time 
behavior of the perturbation. 

3 Bloch Decompositions and Modulational Stability 

In this section, we detail the methods of our modulational stability analysis by utilizing 
Bloch-wave decompositions of the linearized problem. In particular, our goal is to rigorously 
justify the Whitham expansion described in the previous section without reference to the 
Evans function. To this end, recall from Floquet theory that any bounded eigenfunction v 
of d x C[u] must satisfy 

v( x + T) = e^v(x) 
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for some 7 G [— it, 7r]. Defining e = e(«) = it follows that u must be of the form 



v(x) = e £X P{x) 



for some k E [— 7r, 7r] and some T-periodic function P which is an eigenfunction of the 
corresponding operator 



With this motivation, we introduce a one-parameter family of Bloch-operators L e defined 
in (13. ip considered on the real Hilbert space Lp er ([0, T]). By standard results in Floquet 
theory, the spectrum of a given operator L £ is discrete consisting of point eigenvalues and 
satisfy 



and hence the L 2 (R) spectrum of the linearized operator <9 x £[it] can be parameterized by 
the parameter e. As a result, the above decomposition reduces the problem of determining 
the continuous spectrum of the operator d x C[u] to that of determining the discrete spectrum 
of the one-parameter family of operators L e . 

As we are interested in the modulational stability of the solution u, we begin our analysis 
by studying the null-space of the un-perturbed operator Lq = d x C[u] acting on T per ([0,T]). 
We will see that under certain non-degeneracy conditions Lq has a two-dimensional kernel 
with a one-dimeisional Jordan chain. It follows that the origin is a T-periodic eigenvalue 
of Lq with algebraic multiplicity three, and hence for |e| <C 1, considering L e as a small 
perturbation of Lq, there will be three eigenvalues bifurcating from the e = state. To 
determine modulational stability then, we will determine conditions which imply these 
bifurcating eigenvalues are confined to the imaginary axis. 

To begin, we formalize the comments of the previous section concerning the structure of 
the generalized null-space of the unperturbed linearized operator Lq = d x C[u] by recalling 
the following lemma from \BiJ\ IBrJKj . 

Lemma 3.1 ( [Br Jj. lBrJK| ). Suppose that u(x;(Io,Eq,cq) is a T-periodic solution of the 
traveling wave ordinary differential equation ()2. 1 j) . and that the Jacobian determinants Te, 
{T, M} aj E, and {T, M, P} a ,E,c o-re non-zero at (ao, Eq,cq) G T>. Then the functions 



(3.1) 



L £ = 



e- £X d x C[u]e< 




e 



<P0 



{T, u} aj E 



1 



1 




'X 




'2 



{u,T,M} a , E ,c 



'U'2 



{T,M} E ,c + {T,M} atE u 



satisfy 



L (j)Q = L Q (f)i = 
Lq4>2 = ~<t>i 



Lq^I = V>2 
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In particular, if we further assume that {T, M} a ^E and {T, M, P} a ,E,c are non-zero at 
(do, Eq, cq), then the functions {<^}? = i forms a basis for the generalized null-eigenspace 



Throughout the rest of our analysis, we will make the non-degeneracy assumption that 
the quantities Te, {T, M} a ^E, and {T, M, P} a ,E,c are non-zero (see the previous section). 
Lemma 13.11 then implies, in essence, that the elements of the T-periodic kernel of the 
unperturbed operator Lq are given by elements of the tangent space to the (two-dimensional) 
manifold of solutions of fixed period and fixed wavespeed, while the element of the first 
generalized kernel is given by a vector in the tangent space to the (three-dimensional) 
manifold of solutions of fixed period with no restrictions on wavespeed. It immediately 
follows that the origin is a T-periodic eigenvalue (corresponding to e = 0, i.e. k = 0) 
of Lq of algebraic multiplicity three and geometric multiplicity two. Next, we vary k in 
a neighborhood of zero to express the three eigenvalues bifurcating from the origin and 
consider the spectral problem 



for |e| <C 1, where we now make the additional assumption that the three branches of the 
function /i(e) bifurcating from the /j,(0) = state are distinclH. Our first goal is to show that 
the spectrum fi(e) and hence the corresponding eigenfunctions v(e) are sufficiently smooth 
(C 1 ) in e. The stronger result of analyticity was proved in |BrJ| using the Weierstrass- 
Preparation theorem and the Fredholm alternative. Here, we follow the methods from 
[JZ5| to offer an alternative proof which is more suitable for our methods. 

Lemma 3.2. Assuming the quantities {T, M} a ^E and {T, M, P} a ,E,c are non-zero, the 
eigenvalues fij(e) of L E are C 1 functions of e for \e\ <1. 

Proof. To begin notice that since [i = is an isolated eigenvalue of Lq, the associated total 
right eigenprojection R(0) and total left eigenprojection L(0) perturb analytically in both 
e (see [K]). It follows that we may find locally analytic right and left bases {vj(e)}^ =1 and 
{vj (e)}| =1 of the associated total eigenspaces given by the range of the projections R(e) 
and L(e) such that Vj(0) = 4>j and Vj(0) = ijjj. Further defining the vectors V = (v 1,^2,^3) 
and V = (vi,V2,vs)* , where * denotes the matrix adjoint, we may convert the infinite- 
dimensional perturbation problem for the operator L £ to a 3x3 matrix perturbation problem 
for the matrix 



In particular, the eigenvalues of the matrix M £ are coincide precisely with the eigenvalues 
Hj{e) of the operator L £ lying in a neighborhood of fi = 0, and the associated left and right 
eigenfunctions of L e are 




L £ v(e) = fi(e)v(s) 



(3.2) 




fj = Vvjj and fj = WjV* 



2 In the case where two or more branches coincide to leading order, more delicate analysis is required 
than what is presented here. In particular, all asymptotic expansions must be continued to at least the next 
order in order to appropriately track all three branches. 
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where Wj and Wj are the associated right and left eigenvectors of M e , respectively. Thus, to 
demonstrate the desired smoothness of the spectrum of the operator L e near the origin for 
|e| <C 1, we need only demonstrate the three eigenvalues of the matrix M £ have the desired 
smoothness properties, which is a seemingly much easier task. 
Next, we expand the vectors 

(3.3) Vj(e) = <pj + e qj (e) + 0(\e\ 2 ) 

(3.4) vj(e) =i> j + eq j (e) + 0{\s\ 2 ) 

and expand the matrix M e as 

M £ = M + eM 1 + O(\e\ 2 ). 

Notice there is some flexibility in our choice of the functions qj and q~j, a fact that we will 
exploit later. Now, however, our calculation does not depend on a particular choice of the 
expansions in (|3.3j) . From Lemma 13.11 a straight forward computation shows that 

/ \ 

(3.5) M = (i/ti,L <h) 

\ / 

where (^i,L &) = \{T, M} a , E {T, M, P} 

a,E,c is non-zero by assumption reflecting the Jor- 
dan structure of the unperturbed operator L . By standard matrix perturbation theory, 
the spectrum of the matrix M £ is C 1 in e provided the entries [Mi] 12 and [Mi^^ of the 
matrix Mi are both zero. Using Lemma 13.11 again, it indeed follows that 

[Mi]i, 2 = (V'o, ^i0i + L ogi ) + (g 0j Lo^i) = 
[Mi] 3i 2 = {1P2, Ltfa + L qi) + {q 2 ,L (f>i) = 

and hence the spectrum is C 1 (IR) in the parameter e near e = as claimed. □ 

As a result of Lemma 13.21 the associated (non-normalized) eigenfunctions bifurcating 
from the generalized null-space are C 1 in the parameter e in a neighborhood of s = 0. 
Moreover, our overall strategy is now clear: since the eigenvalues of M e correspond to the 
eigenvalues of the Bloch operator L £ near the origin, we need only study the characteristic 
polynomial of the matrix M £ near e = in order to understand the modulational stability 
of the underlying periodic wave. However, notice by equation (|3.5p the unperturbed matrix 
Mo has a a non-trivial Jordan block, and hence the analysis of the bifurcating eigenvalues 
must be handled with care. In order to describe the breaking of the two-by-two Jordan 
block described in Lemma 13. II for e small, we rescale matrix M £ in (I3.2D as 

M £ = e- 1 S(E)- 1 M e S(e) 



where 
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In particular, the matrix M £ is an analytic matrix- valued function of e, and the eigenvalues of 
M £ are given by Vj(e) = e~ 1 Hj(e), where Hj(e) represents the eigenvalues of M E . Notice that 
the second coordinate of the vectors in C n+1 in the perturbation problem (I3.2p corresponds 
to the coefficient of u x to variations if) in displacement. Thus, the above rescaling amounts 
to substituting for ip the variable \e\ip ~ vp x of the Whitham average system. Our next goal 
is to prove that the characteristic polynomial for the rescaled matrix M e agrees with that of 
the linearized dispersion relation A(/i,k) corresponding to the homogenized system (I2.12p . 
i.e. we want to prove that 



det ( sM e - ii ( 5(e) _1 y e *y e 5(e) ) = A(>, k) + 0(\n\ A + \k 



4\ 




for some non-zero constant C. To this end, we begin by determining the required variations 
in the vectors V £ and V £ near e = contribute to leading order in the above equation. We 
begin by studying the structure of the matrix M e . 

Lemma 3.3. The matrix rescaled M £ can be expanded as 



M £ = e~ 1 (S-\e) ( Vi I M<£o, 01 + £<?!, 2 ) 5(e) ) 

in a neighborhood of e = 0. In particular, the only first order variations of the vectors V £ 
and V £ which contribute to leading order are q2 and q\. 

Proof. The idea is to undo the rescaling and find which entries of the unsealed matrix M £ 
contribute to leading order. To begin, we expand the non-rescaled matrix M £ from f|3.2[) as 

M £ = M + eMi + e 2 M 2 + 0(\e\ 3 ) 

and notice that Mq was computed in ()3.5p and was shown to be nilpotent but non-zero, 
possessing a nontrivial associated Jordan chain of height two. Using Lemma 13. II a straight 
forward computation shows the matrix M\ can be expressed as 

(i/) , L^o + L q ) (V>o,£i02) + (qo,L (p2) 

Mi = | * (^i,Li0! + L qi) * 

(V>2,£i0o + L q ) {ip2,Li4> 2 ) + (q2,Lo4> 2 ) 

(lp ,Ll<f) ) (^0,^102) + (90, ^002) 

* (ip 1 ,L 1 (f>x+ L qi) * 

(</>2,£i0o) {ip2, L^) + {q2, Lo4> 2 ) 

where the * terms are not necessary as they contribute to higher order terms in the rescaled 
matrix M £ . Similarly, the relevant entries of the matrix M 2 are given by 

* (^0,^201 + £i<7i) + (go,£i0i + £o<7i) * 
M 2 = I * * * 

* (ip 2 ,L 2 (j) 1 + Liqx) + (q 2 ,L l (f) 1 + L qi) * 
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Therefore, it follows the relevant entries of the matrix M £ up to 0(|e| 2 ) can be evaluated 
using only the variations q\ and q\ , as claimed. □ 

The main point of the above lemma is that in order to compute to projection of the 
operator L £ onto the eigenspace bifurcating null-space at e = to leading order, you only 
need to consider the variations in the bottom of the left and right Jordan chains: all other 
variations contribute to terms of higher order. Our next lemma shows the variation in 
the ifjQ direction is also needed to compute the corresponding projection of the identity to 
leading order. 

Lemma 3.4. Define the matrix I £ := /s(E)~ 1 V*V £ S(e)\. Then I £ can be expanded near 
e = as 

I £ = (S- 1 (e){ fa \{<h,<h. + eqi,<h)S(e)) +o(l). 



ip2 + £Q2 



L 2 pe r([0,T]) 



Proof. As in the proof of Lemma I3.3[ we undo the rescaling and find the terms which 
contribute to leading order. To begin, we expand the matrix I £ := (v*V E \ as 

I £ = I + eh + 0(\e\ 2 ) 
Using the fact that (fa, <f>j) = if i ^ j, it follows that 

(fa,4> ) 

o (fa, fa) o 

{ip 2 ,cp 

* {^o,1i) + (<Zo,</>i) * 



2 



and 



* (1P2, qi) + (<Z2, <f>i) * 

from which the lemma follows by rescaling. □ 

With the above preparations, we can now project the operator L £ — /x(e) onto the total 
eigenspace bifurcating from the origin for small |e| <C 1. Our claim is that projected and 
rescaled matrix agrees with the Whitham system, up to a similarity transformation. To 
show this, we begin by showing the characteristic polynomial of the corresponding rescaled 
matrix agrees with the linearized dispersion relation A(/x, k) corresponding to the linearized 
Whitham system. Once this is established, the fact that the bifurcating eigenvalues are 
distinct will imply the desired similarity by the analysis in [JZlJ. 

In order to compute the characteristic polynomial of the matrix projection of the oper- 
ator L £ — /i(e), we must now make a few specific choices as to the variations in the vectors 
V £ and V £ . As noted in the proof of Lemma E21there is quite a bit of flexibility in our choice 



3 BLOCH DECOMPOSITIONS AND MODULATION AL STABILITY 



15 



of the expansion (|3.3p . As a naive choice, we could require that Vj(e) be an eigenfunction 
of L £ for small e. For example, we may choose vq and vi to satisfy 

LoQj = (Ai - Li)<t>j, j = 0,2 

where Ai is the corresponding eigenvalue bifurcating from the origin. Similarly, we may 
choose vo and v% to satisfy 

Llq j = (Xl-L 1 )iP J , j =0,1 

where * denotes complex conjugation. However, this is not necessary: we only need Vj(e) 
to provide a basis for the null-eigenspace of the operator L E . To illustrate this, instead of 
choosing the variation in the 4>\ direction to satisfy an eigenvalue equation we notice there 
are two distinct eigenvalues with expansions 

Hx(e) = eAi + o(|e|), 
fix(e) = eAi + o(\e\) 

corresponding the the eigenfunctions <p\ + eg\ + o(\e\) and 4>i + eg% + o(|e|). In particular, 
it follows that the functions g± and g\ satisfy the equations 

Logi = (Ai - Li)4>i, 
Loh = (Ai - L\)(j>i. 

Defining q\ = (\\ — {^igi — ^x9i \ then, it follows that the function q\ to satisfies 

Lm = (Ai - Ax)"' (Xi(Ai - Li)^i - Ai(Ai - Li)0i) 
= -Li4>\ 

and hence may be chosen to satisfy q\ _L span{-0o ; V'l}- To find a closed form expression for 
qi, notice that 

L (xu x ) = 2u xxx = -L\u x . 

Moreover, a direct calculation shows that the function ue satisfies C[u]ue = and the 
function 

T 



4> = ~{T, M} a<E \ XU X + jT^e 

is T-periodic. In particular, Lq(J) = —L\<\>\ and <p _L spanl^OiV 7 !} an d hence we may choose 
qi = (j). Note there are many other choices for q\ that are possible: our choice is made to 
simplify the forthcoming calculations. Similarly, since ip2 is at the bottom of the Jordan 
chain of the null-space of Lq, we may choose g\ to satisfy the equation 

A)<72 = -L[tp2- 
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Unlike the variation in (f>\, we do not need a closed form expression for g\\ the above defining 
relation will be sufficient for our purposes. 

With the above choices, all the necessary inner products described in Lemmas 13.31 and 
3.41 may be evaluated explicitly. Indeed, straight forward computations show that 


,1/n-l \{T : M}a tE {T,M,P} atEtl 


T E T * 

A/ 1 - | * * 

T(T E {T,M} E , c + Ta{T,M} a , E ) * 

T{T,K} a , E 
} h = * * 

\ * T{T, M} E<C {T, K} a>E + T{T g KE (T a {T, K} a , E - T{T, M} a>E ) 



* \ 

* 

* 



and 
h = 



{T, M} a)E 

-\{T,M} a)E {T,M,P} a ^ t 










o 

±{T,M} atE {T,M,P} atEtC 
-{T,M} E>c T-{T,M} aiE M 

* 2{T, M} atE {K, T, M} a , E , c - {T, M} 2 Ec T - 2{T, M} B>C {T, M} a , E M - {T, M} 2 aE P 

We can thus explicitly compute the rescaled matrices M £ and I e in terms of the underlying 
solution u, which yields the following theorem. 

Theorem 2. Let (ao, Eq,co) 6 T> and assume the matrices ^^q^'^ and g^ E e ) are 
invertible at (ag, Eq,cq). Then the linearized dispersion relation A(p,, k) in (|2.13p satisfies 

det (eM £ - n ls(e)- l V*V £ S(s 



= CA(/i, K ) + 0(H 4 + M 4 ) 

for some constant 0. That is, up to a constant the linearized dispersion relation for the 
homogenized system (I2.12D accurately describes the low-frequency behavior of the spectrum 
of the Bloch- operator L^. 

Proof. A straightforward computation using the above identities implies 

T 3 {T,M}3 



det [eM E -n{ S(e)-%*V e S(e 



({T,M} Eia + T*-2T c T E ) 



2T E 

+ CA{n, K ) + 0(\fi\ 4 + \ K \ 4 ) 
for some non-zero constant C = C(a,E,c). Moreover, using the identity 

u 2 du 



2T C = P E = M a = 



r (E -y(u;a,c)) 3/2 ' 
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which immediately follow from the integral formulas (|2.3p - (|2.5p . along with the fact that 
T a = M E by IpTf]) . it follows that 

T 2 a - 2T C T E = T a M E - T E M a = {T, M} a>E = —{T, M} Ea 

and hence 

det [eM £ - /i ^(^"^/^(e))) = CA(fi, k) + 0(\n\ 4 + |k| 4 ) 

as claimed. □ 

Theorem [2] provides a rigorous verification of the Whitham modulation equations for 
the gKdV equations. While this has recently been established in [JZlj . the important ob- 
servation here is that the verification was independent of the restrictive Evans function 
techniques. Thus, the above computation can be used as a blue print of how to rigorously 
justify Whitham expansions in more complicated settings where the Evans function frame- 
work is not available. As a consequence of our assumption that the eigenvalues of M E be 
distinct for < \e\ <C 1, it follows that there must exist a similarity transformation between 
the matrix 

d £ (eM £ - n (s{e)- l V:V e S{e))) | £=0 

and the matrix arising from the linearized Whitham system (I2.13D . Thus, the variations 
predicted by Whitham to control the long-wavelength stability of a periodic traveling wave 
solution of the gKdV are indeed the variations needed at the level of the linearized Bloch- 
expansion. 



4 Conclusions 

In this paper we considered the spectral stability of a periodic traveling wave of the gKdV 
equation to long-wavelength perturbations. Recently, this notion of stability has been the 
focus of much work in the context of viscous systems of conservation laws [OZ1L ]OZ3, 
IOZ41 ISell IJZ44 IJZ5] and nonlinear dispersive equations [Br J, Br JKl |J2~1 IJZj . While much 
of this work has utilized the now familiar and powerful Evans function techniques, our 
approach of using a direct Bloch-decomposition of the linearized operator serves to provide 
an elementary method in the one-dimensional setting considered here as well as a robust 
method which applies in the more complicated setting of multi-periodic structures. As such, 
our hope is that the simple and straightforward analysis in this paper will be used as a blue 
print of how to justify the Whitham modulation equations in settings which do not admit 
a readily computable Evans function. In future work, we will report on the application of 
this method to the long-wavelength stability of doubly periodic traveling waves of viscous 
systems of conservation laws. 

As a future direction of investigation, we note that Theorem [2] and the assumption that 
the branches of spectrum bifurcating from the origin are distinct suggests a possibly more 
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algorithmic approach to justifying the Whitham modulation equations via the above Bloch- 
expansion methods which furthermore does not rely on comparing the linearized dispersion 
relations as in Theorem [2] it should be possible to justify the Whitham expansion by direct 
comparison of inner products and showing the mentioned similarity directly. As a first 
step, we suggest beginning this line of investigation by studying the gKdV equation as in 
the current paper and [TZl] and expressing the full Whitham system in terms of inner 
products. Demonstrating the desired similarity in this way, while not completely necessary, 
could substantially simplify the amount of work required to justify the Whitham equations 
by not requiring an initial justification at the spectral level (a possibly daunting task in 
more general situations). 
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